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Abstract 

Possibilities of discriminating neutrons and y rays in the AGATA y-ray tracking spectrometer have been investigated 
with the aim of reducing the background due to inelastic scattering of neutrons in the high-purity germanium crystals. 
This background may become a serious problem especially in experiments with neutron-rich radioactive ion beams. 
Simulations using the Geant4 toolkit and a tracking program based on the forward tracking algorithm were carried 
out by emitting neutrons and y rays from the center of AGATA. Three different methods were developed and tested 
in order to find "fingerprints" of the neutron interaction points in the detectors. In a simulation with simultaneous 
emission of six neutrons with energies in the range 1-5 MeV and ten y rays with energies between 150 and 1450 keV, 
the peak-to-background ratio at a y-ray energy of 1.0 MeV was improved by a factor of 2.4 after neutron rejection 
with a reduction of the photopeak efficiency at 1.0 MeV of only a factor of 1.25. 

Key words: Gamma-ray tracking, AGATA, HPGe detectors, neutron-gamma discrimination, inelastic neutron 
scattering, Geant4 Monte Carlo simulations 
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1. Introduction 

One of the large challenges of contemporary experi- 
mental nuclear structure physics is to study exotic nu- 
clides, which lie far from the line of yS stability. For pro- 
duction of such proton- and neutron-rich nuclei the us- 
age of radioactive ion beams is essential. The intensity 
of such beams is, however, several orders of magnitude 
lower than what is obtainable with stable beams and the 
production rates of the exotic nuclides are very low. In 
order to partly compensate for the low production rates 
highly efficient detection techniques are being devel- 
oped. One such technique is y-ray tracking |1|. Two 
y-ray tracking arrays made of segmented high-purity 
germanium (HPGe) detectors are currently under devel- 
opment, AGATA L2J in Europe and GRETA Li] in the 
USA. These future instruments promise an increase in 
performance over existing detector arrays by several or- 
ders of magnitude for certain experimental conditions. 
This increase in the experimental sensitivity opens up 
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the opportunity to study very weak reaction channels 
and interesting nuclear phenomena. 

One of the most common types of nuclear reac- 
tions used in nuclear structure physics is the heavy-ion 
fusion-evaporation (HIFE) reaction in which two nuclei 
fuse to form a compound nucleus that decays by evapo- 
rating a number of light particles, mainly neutrons, pro- 
tons, and/or a particles. In contrast to the light charged 
particles, neutrons can travel long distances and inter- 
act within the germanium detectors together with the y 
rays emitted by the residual nuclei. Neutrons are de- 
tected indirectly in AGATA either via elastic scattering 
with the Ge nuclei or via nuclear reactions which may 
produce neutrons, charged particles and y rays. Elas- 
tic scattering, "'"Ge(n,n)™'Ge, gives the largest contri- 
bution to the total cross-section in the neutron energy 
range of interest for HIFE reactions, from about 1 to 10 
MeV, while inelastic scattering, "^'Ge(n,n'y)"'"Ge, has 
the second largest cross section in this energy range ID. 
The expected probability of neutrons with energies 1- 
10 MeV to be detected by AGATA is about 50 %, if the 
low-energy threshold of the Ge detector signals is set to 
5 keV (target value for AGATA) |5J. 

In several recent articles (H |5] EJ, experimen- 
tal studies of neutron interactions in segmented HPGe 
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detectors were reported with the aim of reducing the 
neutron-induced background in the y-ray spectra and 
to utiUze segmented HPGe detector arrays to measure 
neutrons. In an earlier Monte Carlo simulation Q of 
the AGATA y-ray spectrometer, the background in the 
y-ray spectra due to the inelastic scattering of neutrons 
was discussed and different possibilities of background 
suppression were investigated without any conclusive 
results. The neutron background is expected to become 
a serious problem in future studies, especially when us- 
ing neutron-rich radioactive ion beams of the coming fa- 
cilities, e.g. NUSTAR/FAIR [8 9|, and SPIRAL2 flOll . 
and in experiments where a low background in the y-ray 
spectra is of importance. 

In this work, a Geant4 IfTTII based Monte Carlo sim- 
ulation of the AGATA spectrometer lfT2l has been per- 
formed in order to investigate the possibility of reducing 
the neutron induced background by y-ray tracking. The 
tracking algorithm of the mgt program lfT3lfT4l has been 
used to identify the neutron interaction points and meth- 
ods for eliminating the y rays originated from inelastic 
neutron scattering have been developed with a special 
emphasis on not reducing the photopeak (- full-energy 
peak) efficiency of the y rays of interest. 

An obvious method of distinguishing between neu- 
trons and y rays is to use the time-of-flight technique. 
This possibility was previously investigated in a simu- 
lation of the AGATA array ||4l. It is still an open ques- 
tion if the time resolution that can be achieved by us- 
ing pulse-shape analysis with segmented HPGe detec- 
tors will be sufficient to discriminate neutrons from y 
rays. The time-of-flight technique may be a comple- 
mentary procedure that can be used together with the 
tracking based methods presented in this work. 

2. Monte Carlo simulations and yray tracking 

The interaction points of neutrons and y rays in 
the HPGe detectors of AGATA were simulated with 
the Geant4 toolkit versions 4.9.0 and 4.9.2 using the 
neutron cross section data libraries G4NDL3.11 and 
G4NDL3.13, respectively. Within the low-energy mod- 
els of neutron scattering used by Geant4 prior to version 
4.9.2, errors related to the angular distribution of scat- 
tered neutrons as well as errors in the energy distribution 
of the scattered neutrons, of the recoiling nuclei, and of 
the y rays produced after the inelastic scattering of neu- 
trons, were reported on the Geant4 problem tracking 
system ifTSl . In version 4.9.2 of Geant4 these bugs have 
been corrected. In our work we used for Geant4 ver- 
sion 4.9.0 a modification of the Geant4 code that was 
suggested by J. Ljungvall et al. |4 | and which corrects 



for the errors in the energy distribution of the germa- 
nium recoils. We noted that our results regarding in- 
elastic scattering of 1 to 5 MeV neutrons on germanium 
obtained by Geant4 version 4.9.0, with the correction 
suggested by J. Ljungvall et al. |4| were similar to the 
ones obtained with Geant4 version 4.9.2. 

The full AGATA aiTay with 180 HPGe crystals was 
used in the Geant4 simulations. The AGATA HPGe de- 
tectors were arranged in a An geometry around the cen- 
tral source position (= target position in an in-beam ex- 
periment). In most of the simulations presented in this 
work only the HPGe crystals (no other material) were 
included. In a few of the simulations, the aluminium 
capsules in which the HPGe crystals are mounted, were 
also included. 

In a typical HIFE reaction the majority of the neu- 
trons are emitted in the energy range 1-5 MeV with a 
high-energy tail that extends up to about 15 MeV. In 
this work we either simulated the emission of mono- 
energetic neutrons or neutrons with a flat energy distri- 
bution from 1 to 5 MeV, which covers the majority of 
the neutrons emitted in a HIFE reaction. 

The recoil energies of the Ge ions are not fully con- 
verted into electron-hole pairs. This gives rise to the 
pulse-height defect (PHD) lfT6l that was included in the 
simulations as suggested in ref. I?). The Ge ionization 
energy Ei is expressed as 

Ei^a-E^, fl = 0.21, Z7= 1.099, (1) 

where Zsr is the Ge recoil energy in the HPGe crystal. 

The output of the Geant4 simulation program that 
contains the energy and three-dimensional position of 
each of the interaction points in AGATA was used as 
input to the mgt y-ray tracking program. At this stage 
mgt has no information whether the interaction is due to 
a y ray or a neutron. Thus, the neutron interactions are 
treated in the same way as the y-ray interactions during 
the tracking process. 

Before applying the tracking with mgt the interaction 
points that are located closer to each other than 5 mm 
were packed together Their energies were summed 
and their new position was calculated as the energy- 
weighted average position. The packed interaction en- 
ergies were smeared by a Gaussian distribution with a 
full-width-at-half-maximum (FWHM) corresponding to 
the typical intrinsic resolution of a large volume HPGe 
detector. The interaction positions were then smeared 
by a Gaussian distribution having a FWHM of 5 mm 
at 100 keV. The smearing parameter had an energy de- 
pendence of (interaction point energy) Interaction 
points that have an energy lower than a threshold energy 



of 5 keV were removed. The selection of this value is 
based on measured ifTTI low-energy thresholds of the 
actual AG ATA HPGe detectors. For accepting an event 
in AGATA the deposited energy in at least one crystal 
must be larger than about 30 to 50 keV. Such a trig- 
ger condition was tested and it was found that it had a 
negUgible effect on the discrimination of events due to 
inelastic scattering of neutrons. Therefore, only the 5 
keV energy threshold was used in this work. 

The mgt code uses a forward tracking method which 
is developed by several groups in Europe flT'TSl and in 
the USA |l9il . The method is based on the clusteriza- 
tion of interaction points belonging to the same initial 
photon. The clusters contain interaction points with a 
limited angular range, which is a result of the forward 
peaked Compton-scattering cross section and the fast 
decrease of the mean free path in Ge as the y ray looses 
its energy. An interaction is allowed to be a member of 
several clusters since at the end of the tracking only the 
cluster with the lowest figure-of-merit (FM) value is ac- 
cepted. The FM calculation is based on a comparison 
of the energy of the scattered y rays calculated by us- 
ing the interaction energies and by applying the Comp- 
ton scattering formula to the positions of the interaction 
points. The probabilities for Compton scattering, pair 
production, and photo absorption (photoelectric effect), 
as well as the probability for the y ray to travel the dis- 
tances between the interaction points in a cluster, are 
also included in the FM. 

As an example, the equation for the FM of a cluster 
with two interaction points, the first one due to a Comp- 
ton scattering and the second one due to photo absorp- 
tion, is given by 



FM = 



Eyfi 



(2) 



Here Ey i(en) is the energy of the scattered y ray after 
the first interaction point 



Ey ^i(en) = Eyfi - Ei, 



(3) 



with Eyfl being the sum of all interaction points (- inci- 
dent energy for a fully absorbed and correctly tracked y 
ray) and E^ the energy deposited at the first interaction 
point. Ey iipos) is also the energy of the scattered y ray 
after the first interaction point, but calculated by using 
the Compton formula 



^^^(POS) = Eyfi 



1 + 



£'y,o(l - COS 6) 



(4) 



where the scattering angle of the y ray, 0, is calculated 
by using the position coordinates of the three interac- 
tion points (target position, first and second interaction 
points in this example) and m^c^ is the rest mass of the 
electron. The probabilities for traveUng the distances ri 
and r2 are given by 



PAiri) = e 



PAiri) = e 



(5) 



where A is the mean free path of y rays in Ge. The prob- 
abilities for Compton scattering and photo absorption 
are given by 



omp 



photo — 



O"Comp(£'y,0) 
O" tot(Eyfl) 

o-photo(£'y,i(en)) 
o-tot(£'y,i(en)) 



(6) 



The differential cross section do-]<j^/dQ., which does not 
have an important effect on our tracking results, is the 
Klein-Nishina Compton scattering cross section for un- 
polarized radiation. 

During the tracking process the FM values of all pos- 
sible permutations of interaction points are evaluated. 
The one with the lowest FM value is chosen as an ac- 
cepted tracked y ray, but only of its FM value is < 1 . It 
should be noted that similar but not identical FM def- 
initions as in eq. [2] are given in refs. ifTSl [T9l . These 
three FM equations are implemented as options in mgt 
and they give similar results (Tf] . 

The properties of the data sets used in this work are 
summarized in table [T] Data sets 1, 4, and 5 were pro- 
duced in order to develop the methods for the discrim- 
ination of neutron and y-ray interaction points. Data 
set 6 contains events produced with a flat distribution of 
neutron energies from 1 to 5 MeV, which well enough 
emulates the energy distribution of neutrons emitted in 
a typical HIFE reaction. The results obtained with data 
set 6 were compared to the ones obtained with data set 
2, which contains events generated by 1 MeV y rays 
emitted from the center of AGATA. The energy 1 MeV 
was chosen to represent a typical average y-ray energy. 
Finally, a more realistic simulation was carried out with 
a flat distribution of neutron energies between 1 and 5 
MeV emitted from the center of AGATA together with 
a y-ray cascade representing a "rotational band" popu- 
lated after a HIFE reaction, data set 7. The results ob- 
tained with data set 7 were compared to the results of 
data set 3, which does not contain any neutron emis- 
sions. 



Table 1 : Energies and multiplicities of y rays and neutrons in the simulated data sets used in this work. The Ge isotope enrichment and the used 
pulse-height defect (PHD) coiTection are also given. The y-ray cascade is a "rotational band" with 10 transitions from Ey = 150 keV to 1450 keV 
with t^Ey = 150 keV. 



Data set number 


Ge enrichment 


y-rays 




neutrons 




PHD 


Ey [keV] 


My 


En [MeV] 


M„ 


1 


^^Ge 


596 


1 









2 


"^'Ge 


1000 


1 









3 


"^'Ge 


cascade 


10 









4 


^^Ge 







1 


1 


no 


5 


^^Ge 







1, 3,5 


1 


yes 


6 


""'Ge 







1-5 


1 


yes 


7 


""'Ge 


cascade 


10 


1-5 


6 


yes 



3. Ge recoil energy distributions 

The n + ""'Ge reaction cross section is more or less 
constant for neutron energies in the range 1-5 MeV and 
it is dominated by elastic scattering |4I- Inelastic scat- 
tering has the second largest cross section, being a factor 
of 8 and 1 .3 smaller than elastic scattering at the neutron 
energies 1 MeV and 5 MeV, respectively. In this energy 
interval the cross sections of other reactions are much 
smaller and have in the context of this work been ne- 
glected. In the elastic scattering process, some of the 
neutron energy is deposited in the detector via the re- 
coiling Ge nuclei. These energies give rise to single-hit 
clusters during the y-ray tracking process and they are 
mostly eliminated from the data since their FM value 
is not good enough and they do not qualify as photo- 
absorption points. Inelastic scattering of neutrons is 
more complicated since in this case the recoiling Ge nu- 
clei usually also emit y rays. These y rays cause an 
unwanted background in the y-ray spectra. 

Simulations of the Ge recoil energy distributions 
were performed by using the most recent version of 
Geant4 (version 4.9.2). Distributions for the different 
Ge isotopes as well as for ""'Ge were investigated. The 
shapes of the recoil energy distributions are determined 
by the fairly isotropic angular distributions and by the 
reaction kinematics. For scattering of 1 MeV neutrons 
on ""'Ge, the recoil energy extends from to 57 keV for 
elastic scattering. The recoil energy distributions due 
to inelastic scattering of 1 to 5 MeV neutrons on ""'Ge 
are shown in fig. [T^. The recoil energies due to inelas- 
tic scattering of 1 MeV neutrons extends from to about 
35 keV for the even-even isotopes ^"Ge, ^^Ge, ^'^Ge, and 
^^Ge and to 53 keV for "Ge. 

Figure [TJ) shows the effect of the PHD correction (eq. 
[T]l on the recoil energies. The maximum recoil energy 



after inelastic neutron scattering with En - I MeV is 
only 16.7 keV after the PHD correction. With a low- 
energy threshold of 5 keV, 52 % of the neutron inter- 
action points for 1 MeV neutrons and 5 % for 5 MeV 
neutrons will not be detected after the PHD correction. 
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Figure 1 : Recoil energies of ""'Ge after inelastic scattering of neutrons 
with energies from 1 to 5 MeV, a) without and b) with PHD correction. 



4. Number of interaction points 

In a simplified example, we have investigated 
whether the y-ray tracking code can discriminate be- 
tween two types of y-rays with the same energies, 
namely 596 keV. The first type of y ray is emitted from 
the center of AG ATA (data set 1), while the second type 
is created within a HPGe detector after inelastic scat- 
tering of 1 MeV neutrons on ^"^Ge (data set 4). This is 
illustrated in the inset of fig. |2] To avoid creating y rays 
of other energies than 596 keV from the inelastic neu- 
tron scattering, the AGATA HPGe detectors were in this 
simulation configured to consists only of ^"^Ge isotopes. 
With 1 MeV neutrons one can only populate the first 
excited 2| state at 596 keV in ^'*Ge (the second excited 



state is at 1204 keV), which then de-excites to the 0| 
ground state by emitting a 596 keV j ray. 

The number of interaction points {npt) in ^"^Ge de- 
termined by the tracking algorithm is shown in fig. |2] 
When the incoming particle is a 1 MeV neutron npt is 
increased on the average by one unit compared to the 
events where the incoming particle is a 596 keV y ray. 
This increase in the npt value is due to the neutron inter- 
action point that is included in the mgt cluster together 
with the y-ray interaction points. 



Gamma, = 596 keV 




"0 2 4 6 8 10 12 14 16 18 20 
Number of interaction points, npt 

Figure 2: Number of interaction points in ^^Ge after tracking. Tlie 
spectra are normalized to iiave tiie same number of counts. 

Simulated and tracked y-ray spectra obtained after in- 
elastic scattering of 1 MeV neutrons on ^"'Ge are shown 
for different npt values in fig. |3] Apart from the peak 
at 596 keV, a bump is visible at higher y-ray energies. 
This bump originates from the mgt clusters where the 
interaction point due to the Ge recoil is included in the 
cluster and its energy is added to 596 keV. As npt in- 
creases, it becomes more likely that one of the interac- 
tion points is due to the scattering of neutrons and the 
bump becomes more prominent. At = 1 MeV the 
width of the recoil energy distribution is 35 keV (see 
section|3]l. The FWHM of the bump in fig. [3]is about 65 
keV for npt - 6, which indicates that there may be two 
or even three neutron interaction points in the clusters, 
apart from the y-ray interaction points. 



5. Neutron-yiscrimination methods 



Based on a comparison of data sets 1 and 4, three 
methods to distinguish between two types of y rays, 
were investigated. 
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Figure 3: Gamma-ray spectra due to inelastic scattering of 1 MeV 
neutrons on ''^Ge (data set 4) obtained after y-ray tracking. In sub- 
figures npt = 2 and «pf = 3, an example with three interactions points 
is illustrated. The first interaction point is due to the recoiling Ge 
nucleus after inelastic scattering of the neutron (open circle). The 
second and third interaction points (filled circles) are due to Compton 
scattering and photo absorption, respectively, of the y ray produced in 
the decay of the excited Ge nucleus. If the neutron interaction point is 
not measured or not included in the mgt cluster, the 596 keV peak in 
the npt = 2 spectrum is incremented. Otherwise, if the recoil energy 
is measured and included in the cluster, the bump in the npt = 3 
spectrum is incremented. 



5. 1. Method 1 : Energy deposited in the first and second 
interaction points 

The first method is based on the low energy depo- 
sition in the detectors by the recoiling Ge nuclei after 
inelastic scattering of neutrons. Distributions of the en- 
ergies of the first (£1) and second interaction points (£2) 
for each mgt cluster are shown in fig. |4] For 1 MeV 
neutrons both the £1 and E2 spectra have large abun- 
dances of counts at low energies, below about 40 keV, 
compared to the spectra of 596 keV y rays emitted from 
the center of AGATA. With a condition (gate) of Ei and 
E2 to be > 40 keV one can discriminate a large number 
of the cases in which the incoming particle is a neu- 
tron compared to a y ray. For the second interaction 
point the shape of the spectrum for incoming 596 keV 



J rays is somewhat different than for the first interac- 
tion point, in particular at low energies. With a gate at 
El > 40 keV slightly more of the photopeak efficiency 
of the 596 ke V peak is lost as compared to the same gate 
on E]. In the simulations presented below the gate on 
£2 was therefore set at a smaller value of the interaction 
point energy. The energies of the third interaction point 
were also studied but they were not used for gating in 
this work since they caused an even bigger loss in the 
photopeak efficiency than E2. 



investigation in a similar way that was used in ref ISj. 
As illustrated in fig. |5] 6g is the geometric angle be- 
tween the line that goes from the center of the detector 
array to the first interaction point and the line that con- 
nects the first interaction point to the second one. Thus, 
6„ is determined from the positions of the interaction 
points relative to the center of AGATA. The angle 6^ is 
the scattering angle calculated from the Compton scat- 
tering formula using the total energy deposited by the y 
ray (sum of aU energies in the cluster) and the energy of 
the first interaction point obtained after y-ray tracking. 
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Figure 4: Energy deposited in a) the first and b) tlie second interaction 
points. Tlie spectra were obtained by sending 596 keV y rays (dashed 
Unes, data set 1) and 1 MeV neutrons (solid lines, data set 4) from the 
center of AGATA. The HPGe detectors were made of ^''Oe isotopes 
in these simulations. The spectra in a) and b) have been normalized to 
have the same number of counts. 




o 1500 
U 



1000 



500 



- E = 


r 0^ 


y 

- 596 keV 


fci: 




npt > 2 




-150-100 -50 50 100 150 



Figure 5: Distribution of the difference AS = - 6^ obtained when 
emitting 596 keV y rays (lower, data set 1) and 1 MeV neutrons (up- 
per, data set 4) from the center of AGATA. The HPGe detectors were 
made of '"^Ge isotopes in these simulations. The angle is obtained 
from the position of the interaction points whereas 6^ is obtained from 
the Compton scattering formula. 



5.2. Method 2: Difference in incoming direction of the 
yay 

The second method is based on the difference in the 
incoming direction of the initial y rays. If they come 
from random positions in space, they must be produced 
by inelastic neutron scattering in the HPGe detectors. If 
they come from the center of the array, they are likely 
not produced by reactions in the detectors. Two differ- 
ent angles, 9g and are extracted and compared in this 



If the y ray comes from the center of AGATA we have 
0g = 0c and we expect to see a sharp distribution of the 
difference A6 = 0g- 9^ centered at 0°. The width of this 
distributions is due to the errors in the determination of 
the energy and position of the interaction points. Such a 
distribution is seen in the lower part of fig. |5] The upper 
part of the same figure shows the distribution of y rays 
produced by inelastic scattering of 1 MeV neutrons on 
^■^Ge. This distribution is wider than the first one and it 



is asymmetric around 0°. The asymmetry is due to the 
low -energy interaction points (due to Ge recoils), which 
are accepted by the tracking code as forward Compton 
scattered j rays with a small angle giving an excess 
of positive A0 values. In this work we have used the 
condition A0 > 15°, in order to discriminate neutrons 
and y rays. If A0 is larger than 15° the cluster is as- 
sumed to contain a y ray from inelastic neutron scatter- 
ing. It should be mentioned that the neutron rejection 
obtained by this method can be improved with a better 
interaction position resolution. An increased position 
resolution leads to a smaller width of the distribution in 
the lower part of fig. |5] which allows the use of a lower 
value of the gate on A0. 
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5.3. Method 3: Selection based on the FMalue 

The third method is based on a comparison of figure- 
of-merit values {FM) of the clusters determined by the 
mgt tracking algorithm when the initial particles are 1 
MeV neutrons and 596 keV y rays emitted from the cen- 
ter of AGATA. In the mgt program a cluster is accepted 
only if all its interactions are valid as one of the three 
y-ray processes: Compton scattering, photo absorption, 
or pair production. The validity depends on the energies 
deposited at each interaction point, the probability of the 
process to happen and the interaction length before the 
process takes place. A cluster that contains a neutron 
interaction point may be accepted {FM < 1) only if 
the neutron interaction point accidentally qualifies as a 
y-ray interaction point. Such a cluster is, however, ex- 
pected to give a higher FM value (eq. [2]i as compared 
to a cluster with only y-ray interactions. 

Figure |6] shows a comparison of the FM values of the 
clusters constructed by the mgt code after the emission 
of 596 keV y rays from the center of AGATA and the 
FM values of the clusters constructed after the emission 
of 1 MeV neutrons. In this figure both "accepted" and 
"good" y rays belong to clusters with FM < 1. The 
energy of the "good" y rays, however, has the value of 
596 keV, i.e. it is equal to the incoming y-ray energy. 
The FM distribution of the "good" y rays is strongly 
peaked at small FM values and it decreases much faster 
than the distribution due to 1 MeV neutrons. With the 
help of fig. |6]we decided to use the gate 0.05 < FM < 1 
in order to discriminate between neutrons and y rays. 

6. Comparison of interaction distances 

We have also investigated and compared the distances 
from the center of AGATA to the first interaction point 
for 596 keV y rays emitted a) from the center of AGATA 



Figure 6: Figure-of-merit (FM) values of the mgt code when the ini- 
tial particles are 596 keV y rays (dashed lines, data set 1) and 1 MeV 
neutrons (solid line, data set 4) emitted from the center of AGATA. 
The HPGe detectors were made of ^^Ge isotopes in these simulation. 
The label "accepted" y rays indicates the y rays which fulfill the re- 
quirement FM < 1, whereas the "good" y rays in addition have an 
energy of 596 keV. The histogram of the 1 MeV neutrons and the 
"accepted" 596 keV y rays have been normalized to have the same 
number of counts. 



(data set 1) and b) from inelastic scattering of 1 MeV 
neutrons (data set 4) on ^'^Ge. The comparison was done 
both for the real first interaction point as obtained from 
the Geant4 simulation (before tracking) and for the first 
interaction point as obtained by mgt (after tracking). 
The results are shown in fig. [7] For 596 keV y rays 
emitted from the center of AGATA, the distributions of 
the distances to the first interaction point before and af- 
ter tracking are very similar. The centroid of the distri- 
bution is at 26 cm, which corresponds to the mean free 
path of 596 keV y rays in Ge, 2.5 cm, plus the 23.5 cm 
distance from the center to the front face of the Ge de- 
tectors. For 1 MeV neutrons there is an abundance of 
events for which the real first interaction distances (be- 
fore tracking) are well above 26 cm (see fig. |7]i. How- 
ever, the events with too large interaction distances are 
eliminated by mgt because the P,i{r\) value becomes too 
small, which leads to an FM value larger than 1 (see 



eqs. 

We have found that setting dififerent conditions on the 
interaction distances in mgt rejects rather few neutron 
events and leads to a relatively large loss of photopeak 
efficiency, compared to what we have achieved with the 
other three methods mentioned above. Since this is also 
true for higher neutron energies (up to 5 MeV), it is con- 
cluded here that gating on the interaction distances is 
not beneficial for improving the neutron-y discrimina- 




Figure 7: Interaction distances measured from tlie center of AGATA 
to the first interaction point for 596 keV y rays (data set 1) and for 
inelastic scattering of 1 MeV neutrons (data set 4) on ^''Ge. Tlie liis- 
tograms have been arbitrarily normalized. 



7. Applications of the neutron-yiscrimination meth- 
ods 

The results of the apphcations of the different 
neutron-y discrimination methods, as described in sec- 
tion [5] are presented in this section. 

7.1. 596 keV yays and 1 MeV neutrons (data set 1 and 
4) 

Fig. [8] shows tracked y-ray spectra in the region 
around the 596 keV peak and its bump obtained from 
events generated by inelastic scattering of 1 MeV neu- 
trons on ^"^Ge (data set 4). The dashed and dot- 
ted lines indicate what is left after applying different 
gates. Quantitative values of the achieved reductions 
are shown in table |2] together with the reduction of the 
photopeak efficiency of the 596 keV peak for 596 keV 
y rays emitted from the center of AGATA (data set 1). 
The aim of the gating is to reduce the 596 keV peak 
and its bump in the neutron induced spectra as much as 
possible while keeping the photopeak efficiency of the 
596 keV peak in the y-ray induced spectra as large as 
possible. 

By using gate set II (£1, E2, A6, FM) the total spec- 
trum is reduced by 69 %, the 596 keV transition by 47 
% and the bump by 92 % in the 1 MeV neutron induced 
spectrum. This causes a loss of 22 % in the 596 keV 
full-energy peak of the spectrum obtained by sending 



Figure 8: Reduction of the 596 keV peak and its bump after apply- 
ing different neutron-y discrimination methods. The spectra were ob- 
tained by sending 1 MeV neutrons on ''^Ge (data set 4). See caption 
of table|2]for the definition of the gates I and II. 

596 keV y rays form the center of AGATA. The two 
neutron-y discrimination methods that contribute least 
to the reduction of the photopeak efficiency of y rays 
emitted from the center of AGATA are the Ey and FM 
tests. If we use an or of gates based on these two meth- 
ods (gate set I) the total neutron induced spectrum is 
reduced by 58 % while the peak is reduced by 32 % and 
the bump by 77 %. This causes a loss of 1 1 % of the 596 
keV full-energy peak in the y-ray induced spectrum. 

Note that the three methods for neutron-y discrimina- 
tion do not exclude each other and there is an overlap of 
neutron events rejected by the different gates. For ex- 
ample, the application of gate set I gives a total neutron 
rejection of 58 %, which is lower than the sum of the 
results given by the E\ and FM gates separately. It is a 
different situation when the same gate set is applied on 
the 596 keV y rays emitted from the center of AGATA. 
In this case there is no overlap of the events rejected by 
the different methods and the total loss in the full-energy 
peak is a sum of the losses obtained by the £1 and FM 
gates separately, as is seen in the last column of table |2] 

7.2. Pulse-height defect correction 

The energies of the neutron interaction points are 
reduced considerably when the PHD correction is in- 
cluded in the simulations, see section |2] If the interac- 
tion energy is below the low-energy threshold of 5 keV 
it will not be measured and the rejection methods, in 
particular method 1 which uses E\ (and Et), will not 
give as good results as the ones shown in table [2] 



Table 2: Reduction of counts in percent of tlie 596 keV peak, its bump, and of the total spectrum after applying different neutron-y discrimination 
gates. The results were obtained by sending 1 MeV neutrons (data set 4) and 596 keV y rays (data set 1) on '^Ge. The gates used were E\ < 40 



keV, £2 < 30 keV, Ae > 15° 


0.05 < FM < 1. The gate combinations (£1 , FM) and (£1 , 


El, Ae, FM) 


are labelled I and II, respectively, and were 


created by using a logical or of the individual gates. 










1 MeV neutrons (data set 4) 




596 keV y rays (data set 1) 


Gates 


596 keV peak 


Bump 


Total 


Peak 


£1 


6 


63 


40 


6 




26 


65 


46 


11 


FM 


28 


55 


44 


5 


I: (£1, FM) 


32 


77 


58 


11 


II: (£1, £2, Ae, FM) 


47 


92 


69 


22 



In figs. |9^) and b) the effect of the PHD correction 
on the y-ray spectra obtained by inelastic scattering of 
1 MeV neutrons on ^'^Ge is shown (data sets 4 and 5). 
Apart from the reduced energies one can also observe a 
reduction of the bump and a corresponding increase of 
the counts in the 596 keV full-energy peak. This is due 
to the fact that some counts are moved from the bump to 
the peak if the recoil energy is not measured. As shown 
in figs. ^), c) and d), the bump gets wider and moves 
to higher energies as the neutron energy increases from 
1 MeV to 5 MeV (cf. the recoil energy distributions 
shown in fig. [T]). When the bump moves to higher ener- 
gies less of the counts in it will be lost due to the 5 keV 
energy threshold. 
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Figure 9: Tracked y-ray spectra obtained after inelastic scattering of 
neutrons on ^''Ge (data sets 4 and 5). The 596 keV and 608 keV peaks 
ai'e due to the 2| — > 0| and 2J — > 2* transitions in ^^Ge, respectively. 



The PHD correction was applied in all simulations 
presented in the rest of this section. 

7. 3. 1 MeV yays and 1-5 MeV neutrons ( data sets 2 and 
6) 

In this subsection the results of inelastic scattering 
of 1-5 MeV neutrons on ""'Ge are presented. The neu- 
trons had a flat energy distribution and the PHD correc- 
tion was applied (data set 6). The different neutron-y 
discrimination methods described in subsection |5] were 
tested and their effect on the photopeak efficiency was 
checked by a simulation of 1 MeV y rays emitted from 
the center of AG ATA (data set 2). 

In fig. 10 the histogram with the solid line shows 
an ungated y-ray spectrum (no neutron rejection gates 
were applied) whereas the dotted histogram shows what 
is left after the application of gate set II (Ei, Ei, AO, 
FM). Quantitative results of the application of the dif- 
ferent gates are shown in table [3] 

With gate set II the total counts in the neutron in- 
duced y-ray spectrum was reduced by 54 %, the 834 
keV peak by 39 %, and its bump by 76 % (see table[3]l. 
This causes a loss of 14 % of the 1 MeV peak in the 
spectrum obtained by sending 1 MeV y rays from the 
center of AGATA (data set 2). If instead, gate set I is 
applied, the total neutron induced spectrum is reduced 
by 40 %, the 834 keV peak by 23 %, and its bump by 
58 % whereas the loss in the 1 MeV full-energy peak is 
only 5 % in the 1 MeV y-ray induced spectrum. 

7.4. Flat energy distribution of 1 to 5 MeV neutrons and 
a yray cascade ( data sets 3 and 7) 

In this subsection the results of applying the neutron 
rejection methods on events in which neutrons are emit- 
ted in coincidence with y rays are presented. A simula- 
tion of simultaneous emission of 6 neutrons and a y-ray 
cascade of multiplicity 10 (data set 7) was compared 



Table 3: Reduction of counts in percent of the 834 keV peak (2]^ — > 0| transition in '^Ge), its bump, and of the total spectrum after applying 
different neutron-y discrimination gates. The results were obtained by sending neutrons with a flat energy distribution from 1 to 5 MeV (data set 6) 



and I MeV y rays (data set 2) 


on ""'Ge. The gates used were Ft 


< 20 keV, £2 < 15 keV, Ae > 15°, 


0.05 < FM < 1. The 


jate combinations (£1, 


FM) and {Ei, Et, A8, FM) are labelled I and II, respectively, anc 


were created by using 


a logical or of the individual gates. 






1 -5 MeV neutrons (data set 6) 




1 MeV y rays (data set 2) 


Gates 


834 keV peak 


Bump 


Total 




Peak 


El 


1 


38 


11 




1 


Ae 


25 


51 


36 




11 


FM 


23 


45 


36 




4 


I; (El, FM) 


23 


58 


40 




5 


II: (£1, E2, A0, FM) 


39 


76 


54 
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Figure 10: Tracked y-ray spectra obtained after inelastic scattering of 
1-5 MeV neutrons on ""'Ge (data set 6). The solid histogram shows 
the y-ray spectrum obtained without any neutron rejection gates while 
the dotted histogram shows what is left after applying gate set II (£1, 
£2, Ae, FM). 



to one with no neutron emission (data set 3). The alu- 
minium capsules of the HPGe detectors were included 
in this simulation. Fig. [TT^ shows the large and com- 
plex background caused by the inelastic scattering of 
neutrons compared to the spectrum for which no neu- 
trons were emitted. 

The influence of the neutrons can be quantified by the 
peak-to-background ratio (PTB) and the total photopeak 
efficiency (e) of the y-ray peaks in the spectra. The re- 
sults are summarized in tables |4] and |5] The PTB ratio 
for a y-ray peak was determined as the number of counts 
in the peak area divided by the background area, which 
was defined as a region of ±2cr around the centroid of 
the y-ray peak, with cr being the width of the peak. The 
two gate sets, I and II, were are used for the neutron- 
y discrimination. The ratios of the PTB values, Rptb, 
which show the improvement of the PTB values when 
applying the neutron rejection gates I and II on the data 
set with M„ = 6, are also given in the table. The ratios 
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Figure 1 1 : a) Tracked spectrum for a y-ray cascade of multiplicity 10 
when no neutrons are present in the events (gray, data set 3) and when 
6 neutrons are emitted in coincidence with the y rays (black, data set 
7). b) Tracked energy of the same y-ray cascade of multiplicity 10 
emitted in coincidence with 6 neutrons before (black) and after (gray) 
neutron rejection. 



RpTB are also plotted as a function of the y-ray energy 



in fig. 12 1. 

As seen in figs. [TT^ and b, the neutron emission as 
well as the neutron rejection method has a major influ- 
ence on the PTB ratios. On the average, the emission 
of 6 neutrons reduces the PTB ratio by a factor of 2.6, 
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Table 4: Peak-to-background ratios (PTB) for the emission of and 6 neutrons in coincidence with a y-ray cascade of multiplicity 10 (data sets 2 
and 7). The RpxB values are ratios of the PTB ratios extracted from the gated and ungated spectra with M„ = 6 using gate set I (column 4/column 
3) and II (column 5/column 3). See caption of table[3|for a description of gates I and II. 

Peak-to-background ratio, PTB ^ptb 



7? 

Ey 


No gates 


No gates 


Gates I 


Gates II 


I 


II 


M„ =0 


M„ = 6 


M„ = 6 


M„ = 6 


100 


32.0 


25.2 


23.3 


23.6 


0.93 


0.94 


250 


16.6 


10.9 


12.0 


12.8 


1.10 


1.17 


400 


34.6 


15.5 


24.3 


26.5 


1.57 


1.71 


550 


36.5 


15.3 


27.6 


30.3 


1.80 


1.98 


700 


35.7 


15.4 


29.7 


35.9 


1.93 


2.33 


850 


35.8 


8.3 


17.9 


22.1 


2.15 


2.65 


1000 


37.5 


13.7 


26.9 


32.8 


1.97 


2.40 


1150 


35.4 


13.7 


26.0 


34.3 


1.89 


2.50 


1300 


57.2 


18.7 


40.1 


49.5 


2.00 


2.46 


1450 


110.2 


25.0 


56.9 


73.2 


2.28 


2.93 
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tron rejection gates I and II, respectively. The 7?^ values 
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Figure 12: Ratio of a) peak-to-background values RpxB and b) photo- 
peak efficiencies i?f with and without neutron rejection for = 6. 
Open and filled circles were obtained by using gate set I and II, re- 
spectively. The data point for the 850 keV y ray was excluded from 
the plots because it is located on the bump of the 834 keV transition 
(see fig. |10) and therefore difficult to analyze. 



plotted in fig. 12 5 are the ratios of the e values with and 
without neutron rejection for gates I and II and M,, = 6. 

By using gate set II, the photopeak efficiency of the 
y-ray peak at 1 MeV was reduced by a factor of 1.25 
(20 %) while the e value of the 1039 keV peak, due to 
the 2| 0| transition in ^"^Ge, was reduced by a factor 
of 2.02 (51 %). With gate set I the e value of the y-ray 
peak at 1 MeV was reduced by a factor of 1.13 (11 %) 
and the 1039 keV by a factor of 1.57 (36 %). 

A test of the dependency of the /?ptb and ratios 
on the low-energy threshold was also performed. At a 
y-ray energy of 1 MeV and with gates II the /?ptb was 
reduced from 2.4 with an energy threshold of 5 keV to 
1.9 with a threshold of 30 keV. The value was not 
effected by the energy threshold. 



whereas the neutron rejection method increases it by a 
factor 1.8 and 2.1 for gates I and II, respectively. At a 
y-ray energy of 1 MeV, the PTB ratio is improved from 
13.7 to 32.8, i.e. by a factor of 2.4 for gates II. The PTB 
ratio at 1 MeV is 37.5 for a simulation with no emitted 
neutrons (Mn - 0). 

The influence of the neutron emission and of the neu- 
tron rejection methods on the photopeak efficiency was 
also investigated. The emission of neutrons have a neg- 
ligible influence on the photopeak efficiency. This can 
be seen by comparing the e values for Mn = 6 to the 
once obtained for Mn = in table [5] The neutron rejec- 
tion, however, decreases the e values. On the average, 
e is reduced by a factor of 1.16 and 1.24 for the neu- 



8. Summary and conclusions 

Monte Carlo simulations of neutrons interacting with 
the HPGe detectors of the AGATA spectrometer were 
performed. The possibility of using y-ray tracking for 
discrimination of interactions due to inelastic scatter- 
ing of neutrons on Ge was investigated with the aim to 
suppress the neutron induced background in the y-ray 
spectra. Three methods were developed to find "finger- 
prints" of the neutron interaction points in the detectors 
and to use them for neutron rejection. The first method 
makes use of the low Ge recoil energies deposited in the 
detectors by the scattered neutrons. The second method 
utilizes the random incoming direction of the initial y 



Table 5: Photopeak efficiency (6) for ffie emission of and 6 neutrons in coincidence wiffi a y-ray cascade of multiplicity 10 (data sets 2 and 7). 
The i?f values are ratios of the photopeak efficiencies extracted from the gated and ungated spectra with M„ = 6 using gate set I (column 4/column 
3) and II (column 5/column 3). See caption of table[3|for a description of gates I and II. 

Photopeak efficiency, e 7?^ 



7? 

Ey 


No gates 


No gates 


Gates I 


Gates II 


I 


II 


M„ =0 


M„ = 6 


M„ = 6 


M„ = 6 


100 


0.60 


0.57 


0.56 


0.56 


0.98 


0.98 


250 


0.52 


0.49 


0.35 


0.34 


0.71 


0.69 


400 


0.48 


0.46 


0.37 


0.34 


0.80 


0.74 


550 


0.44 


0.41 


0.35 


0.32 


0.85 


0.78 


700 


0.41 


0.39 


0.33 


0.31 


0.85 


0.80 


850 


0.38 


0.36 


0.32 


0.29 


0.89 


0.81 


1000 


0.37 


0.35 


0.31 


0.28 


0.89 


0.80 


1150 


0.35 


0.32 


0.29 


0.27 


0.91 


0.84 


1300 


0.23 


0.31 


0.28 


0.26 


0.90 


0.84 


1450 


0.31 


0.30 


0.27 


0.25 


0.90 


0.83 



rays that are produced by inelastic neutron scattering. 
The third method makes use of differences in the figure- 
of-merit values evaluated by the y-ray tracking program 
for y-ray and neutron induced interactions. These meth- 
ods were first tested on a simplified case in which neu- 
trons with energies 1 to 5 MeV and multiplicity 1 were 
emitted from the center of the AGATA spectrometer 
With a combination of the three neutron-y discrimina- 
tion methods, the 834 keV transition in ^^Ge was re- 
duced by 39 % whereas the bump structure that belongs 
to this transition was reduced more effectively, by 76 
%. The bump originates from the tracking of clusters 
where the Ge recoil energy is added to the energy of 
the y ray produced by inelastic scattering of the neu- 
tron. The neutron rejection methods were also tested on 
a case where 6 neutrons with energies between 1 to 5 
MeV were emitted in coincidence with a cascade of 10 
y rays. The influence of the neutron emission and the 
neutron rejection methods were quantified by determin- 
ing the peak-to-background ratios and total photopeak 
efficiencies of the y-ray peaks in the tracked spectra. As 
a result, two sets of gates, which are combinations of 
the three neutron rejection methods, were suggested and 
used. At a y-ray energy of 1 MeV, the first gate set in- 
creases the peak-to-background value by a factor of 2.0, 
whereas the second set of gates increases it by a factor 
of 2.4. The same gates reduce the photopeak efficiency 
at 1 MeV by a factor of 1.13 and 1.25, respectively. The 
future plan is to test the neutron-y discrimination meth- 
ods on real data obtained by the AGATA spectrometer. 
Neutron-y discrimination is also important for measure- 



ments of high-energy y rays e.g. from the decay of giant 
resonances at high temperatures [20]. For this purpose 
it is necessary to investigate neutron interactions in the 
energy interval between 5 to 20 MeV. This work is on- 
going. 
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